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We present a thermodynamic formalism to study the full counting statistics (FCS) of charge 
transport through a quantum dot coupled to two leads in the resonant-level model. We show that 
a close analogue of equilibrium phase transitions exists for the statistics of transferred charge; by 
tuning an appropriate 'counting field', crossovers to different dynamical phases are possible. Our 
description reveals a mapping between the FCS of a given device and current measurements over a 
range of devices with different dot-lead coupling strengths. Further insight into features in the FCS 
is found by studying the occupation of the dot conditioned on the transported charge between the 
leads. 
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Introduction — The last two decades have seen signif- 
icant interest in calculating and understanding the FCS 
of charge transport in mesoscopic systems [THE] ■ Beyond 
the more readily accessible mean and variance, studying 
the FCS of transferred particles gives insights into coher- 
ences [5], particle interactions [TU] and charge fractional- 
ization [TT] . Various theoretical methods exist to extract 
FCS in mesoscopic systems: seminal works introduced an 
auxiliary system which is coupled to the mesoscopic sys- 
tem of interest during quantum transport [I] . An alter- 
native theoretical approach is a two-point measurement 
scheme where the transferred charge is determined by 
the charge difference in a reservoir between two measure- 
ments separated in time [HI EH E3] ■ This latter approach 
gives manifestly positive probabilities for the transferred 
charge and is minimally invasive since the coherent evo- 
lution of the system is not interrupted or altered between 
the two measurements. Over the last few years, experi- 
ments have provided increasingly detailed measurements 
of FCS. For example, a quantum point contact has been 
used to measure the sequential tunnelling of electrons 
through a quantum dot in real time [5J [2] . 

Recent work has developed an "s-ensemble" ap- 
proach |15j for studying the counting statistics of open 
quantum systems |16| which is based on the application 
of large-deviation methods [17] . This approach allows an 
analogy to be developed between quantum trajectories 
and traditional thermodynamics. This large-deviation 
method has been applied to the time record of quantum 
jumps obtained by unravelling a Markovian master equa- 
tion (MME) and has revealed new features in the dynam- 
ics of open quantum systems [18]. Most notable is the 
existence of dynamical phase transitions between phases 
with different quantum-jump rates. These transitions 
may be crossed by tuning system parameters or by ad- 
justing a counting field which biases the system towards 
rare trajectories where the number of quantum jumps dif- 
fers from the mean. This method has also been applied 
to transport problems approximated by a MME [19j . 

In this Letter, we introduce a new formulation of the 
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FIG. 1. A schematic diagram of the quantum dot and semi- 
infinite conducting leads 1 and 2 with thermal occupations 
Hp and at temperature T and chemical potentials /xi 
and fj,2 in the non-equilibrium steady state. 



s-ensemble which captures features beyond those acces- 
sible with a MME description of the dynamics. We use 
the two-point measurement approach to study quantum 
transport, where we consider a trajectory to be speci- 
fied only by the charge transferred between two measure- 
ments separated by a time t. This allows us to study the 
statistics of transport without making Born and Markov 
approximations. We apply the method to an exactly 
solvable model for a quantum dot coupled to two free- 
electron leads and find a rich dynamical phase diagram 
not captured in the MME approach [TS]. The approach 
uncovers a mapping relating the FCS of a system to the 
average current in a class of related systems. This re- 
markable property of the model is not revealed in other 
approaches to FCS. 

Transport model — We study a single-level quan- 
tum dot coupled to two semi-infinite conductors held 
at different chemical potentials, shown schematically in 
Fig. [I] The universal low-energy behavior of free elec- 
trons scattering from a single dot level is captured by 
the resonant-level model p~3|. I20H22] . Using this model 
and the two-point measurement method, Bernard and 
Doyon 13 recently provided an exact derivation of the 
characteristic function for charges transferred through 
the dot after long times, and demonstrated its equiva- 
lence to the established result obtained by Levitov and 
coworkers pQ. The steady-state distribution of charge 
transferred through the dot Aq during a time t takes 
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a large deviation (LD) form P t (Aq) = e-^^il^) [T7]. 
Therefore, the moment generating function Z t (s), found 
via Laplace transform, also takes a LD form Z t (s) = 
£ A? e- sA9 P t (Ag) = e* e « 

This mathematical framework allows a thermodynamic 
analogy to be made where the convex LD functions 
(j)(Aq/t) and 0(s) are considered akin to, respectively, en- 
tropy and free energy densities [HI [16]. Z t (s) plays the 
role of a dynamical partition function, while the param- 
eter s can be considered a izme-intensive conjugate field 
to the izme-extensive number of transferred charges Aq. 
Tuning this 's-field' allows us to bias the system towards 
rare occurrences where Aq is greater (s < 0) or smaller 
(s > 0) than the average. We identify dynamical phase 
transitions with singular regions in 9(s) and establish an 
order parameter, the average current I(s) = (Aq) /t, as 
a way to distinguish dynamical phases. Through its s- 
derivatives, I(s) contains the full set of cumulants for 
charge transport and hence provides a description of the 
FCS of the system. 

After linearising the spectrum about the Fermi en- 
ergy and unfolding the model, we consider the Hamil- 
tonian 



thus given by 
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We consider the case where the system has a steady- 
state density matrix p ss |13l 123) with different chemi- 
cal potentials pi in leads i = 1 and 2 (see Fig. [I]). We 
find the transferred charge by projecting on to charge 
states of lead 1 at times separated by t: after pro- 
jecting on to states of lead 1 with charges qo and q 
at times and t, the density matrix takes the form 



Pile = P q U t P qoPss P qo UjP q , with U t 



lHt the evo- 



lution operator and P q the charge projection operator on 
lead 1. The probability of measuring charge qo in the 
first measurement and then charge q after time t is thus 
Tr(p qqa ) and the LD function for the transferred charge, 
Aq = q - q, is 



8( S )=t- 1 \ogJ2e- s(q "- q) MP qq «)- 
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From this, we extract the order parameter, the aver- 
age current I(s) = —9'(s), where the prime denotes s- 
differentiation. We also study the charge fluctuations 
A/ 2 (s) = [(Aq 2 ) - (Aq) 2 }/t = 6"(s). 
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(V>j(0)d + dty^O^+erfte. Following [T3], 6{s), from its definition in Eq. @, is 
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Here, d^ (d) is the fermionic creation (annihilation) op- 
erator for the dot at energy e and ipj(x) (ipj(x)), with 
j = 1,2 the creation (annihilation) operators for the 
unfolded leads 1 and 2 at position x. The leads have 
length L, which we will take to infinity, and couple to 
the dot with tunnelling matrix element r. Following con- 
vention [22] , we have set the energy scale of the problem 
by choosing unit Fermi velocity, with density of states at 
the Fermi energy g = l/2n. 

In the limit L — > oo, the Hamiltonian (IT]) is diag- 
onalised with even (odd) mode operators aj^(b p ) with 
eigenenergy p such that H = J p(a^a p + bj ) b p )dp. The 
field operators are then expressed 



*(«) = / glog{l + T p [r,W(4 2 )-l)(l-e s ) 

+ f(f -l)(l-e- s )]}, (5) 

where np 1 ^ (n p 2 ^) is the thermal mode occupation at en- 
ergy p for lead 1(2), shown in Fig. [JJ This result matches 
the Levitov formula [TJ , as proven in [T3] . The s-ensemble 
approach also provides a convenient way of investigat- 
ing the state of the dot conditioned on the transported 
charge by finding the dot occupation for rare charge tra- 
jectories. We consider the dot occupation n^(s) in the 
s-biased distribution of measured charges, 
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in the Heisenberg representation. The even mode hy- 
bridises with the dot level such that w p = r 2 /[— t 2 +i(p— 
e)] and a phase shift is introduced into the even mode due 
to scattering off the dot, such that e p = e l5p = 1 — w p 
when x > and e p = 1 for x < 0. The transmission 
coefficient through the dot for electrons with energy p is 



which can be evaluated analytically for the resonant-level 
model (see Supplemental Material [2"1|). 

Dynamical Phases — We now discuss the statistics 
of transferred charge in the s-ensemble picture. To be 
concrete, we fix pi = —30 below the dot energy and 
study the charge dynamics with different s-fields as [i\ 
is changed. A summary of our results is shown in Fig. [2] 
where we set r = 2 and, without loss of generality, choose 
the dot energy e = 0. 

We first examine the zero-temperature dynamical 
phase diagram. Using the mean current I(s) = —9'(s) as 
an order parameter to distinguish different regimes, we 
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FIG. 2. (Color online.) Plots of (a,b) the current order pa- 
rameter I(s), (c,d) current fluctuations AI 2 (s) 1 and (e,f) the 
dot occupation probability n<j(s) for T — (a,c,e) and T = 1.5 
(b,d,f) as a function of s and fi±. We set fj,2 = —30 and e = 
as labelled on the ordinate axis. 

find three distinct 'phases' separated by sharp crossovers 
in both s and Figure [2^a) shows distinct regions 
where charge flows from left to right (I(s) > 0), no 
charge flows (I(s) ~ 0) and charge flows from right to left 
(J(s) < 0). Analogous to equilibrium phase transitions, 
Fig |2](c) shows that crossover lines between dynamical 
phases are marked by sharp peaks in current fluctuations 
A7 2 (s), except where these lines are parallel to the s-axis. 
The phases are divided by distinct crossovers which be- 
come sharper with increasing fii — /i 2 . However, when 
Mi = A*2, n ot only is the average current zero, there are 
no fluctuations in transferred charge. Therefore, appli- 
cation of the s-field is unable to modify the distribution 
of charge trajectories and J(s) = everywhere. Upon 
raising fj,± above fi 2 , charges can flow from the left lead 
to the right. We now see that applying the s-field al- 
lows new dynamical phases to be reached by biasing the 
mean current to larger values if s < and suppressing 
the current for s > 0. 

The dynamical phase structure is modified at finite 
temperature as there is always a finite probability of 
charge flowing in either direction. Figure p[b) illustrates 



the case of T = 1.5. When \i\ = /i2, the system is in equi- 
librium with 1(0) = 0. However, at finite T, there are 
fluctuations in the charge transported between measure- 
ments and application of the s-field takes I(s) through 
sharp crossovers even when /ii — ^2- As /ii is increased, 
dynamical phases with positive and negative 7(s) are still 
present, but with increasing forward bias the s-field re- 
quired to reach the negative current regime grows in pro- 
portion to Hi/T because transport against the bias be- 
comes increasingly rare. The phase crossovers are again 
marked by peaks in the fluctuations A/ 2 (s), but we now 
see finite fluctuations within the phases with with large 
\I(s)\ (see Fig. ^d)). Within these phases, AI 2 (s) is 
virtually independent of s and the potential bias and is 
proportional to T. There are smaller peaks in A/ 2 (s) in 
the low-bias regime /ii < e which are additional to those 
at the phase boundaries. We return to these features 
later. 

The dynamics of the system are constrained by a fluc- 
tuation theorem (FT) for the rate of entropy produc- 
tion [5], and this imposes a structure on the dynam- 
ical phase diagram. The probabilites of transporting 
charge in the forward and reverse directions are related 
by P t (Aq)/P t (-Aq) = e 2s '=°. From Eq. (§, s J=0 = 
(pi — \i%)l1T is the s-field required to bring the current 
to zero. As a consequence of the FT, Bis) is symmetric 
about s/ = o such that I(s — sj = q) = — 7(sj = q — s). The 
form of s/ = o illustrates that the s-field shares character- 
istics with a real potential bias: at finite T it can be 
tuned to oppose the bias /ii — [ii and bring the system 
to equilibrium. 

FCS from current measurements — A remarkable fea- 
ture of the dynamical phases at T = is that rare charge 
trajectories, where s ^ 0, capture the typical charge tra- 
jectories for all dot-lead coupling stengths r. Specifically, 
from Eq. |5]), 

hence there exists a simple mapping between the behav- 
ior at a given s and that of the unbiased (s = 0) dynam- 
ics by modifying the coupling r — > t s = e~ s / 4 r. The 
phase diagram for another coupling r' is thus straight- 
forwardly related to the phase diagram at r by a trans- 
lation along the s-axis where s — > s — 4^1og(r'/r), with 
£ = sgn(/ii — 1x2). Furthermore, since I(s) encodes all the 
cumulants of the current through its derivatives, the FCS 
for a particular r can be inferred from the measurements 
of the average current for a range of coupling strengths. 

More generally, at finite T, we find that FCS of 
a system with an s-bias applied are replicated in the 
s = FCS for another quantum-dot system. In gen- 
eral, energy-dependent couplings must be introduced to 
the model (see Supplemental Material [H]). However, 
provided \^i^\/T 3> 1, with [12 < e < fi\, modifica- 
tions to the form of the Hamiltonian (fll) are not needed. 
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FIG. 3. (Color online.) Plots of I(s) for m = -/t 2 = 30. (a) 
shows the translational property of I(s) with r, at zero and 
finite T, as labelled, (b) demonstrates I(s) (solid lines) can 
be found from the unbiased current (dashed lines) through 
quantum dots with the modified r — > r a and /ti^ — > Mi, 2 ( see 
main text), for r = 1 and different T, as labelled. 



The translation property of I(s) still holds, as shown in 
Fig.[3Ja). Furthermore, after introducing the same mod- 
ified coupling r — > t s and mapping /z^a — > Mi 2> where 
/if 2 satisfies 



Mfo = Mi,2±Tlog 



4(Mi, 



(8) 



we find that we can determine J(s) from values of the 
unbiased current, 1(0), in resonant-level systems with 
coupling t s and chemical potentials /if 2 ■ We demon- 
strate in Fig. [3jb) that good agreement exists between 
this scheme and the exact I(s) for a range of tempera- 
tures when s < si = q. I(s > s/ = o) can be inferred using 
the reflection antisymmetry in I(s) described by the FT. 
This simple structure for the mapping no longer holds 
when the bias is very small as additional features in the 
FCS emerge. These manifest as additional peaks in the 
fluctuations A/ 2 (s) (see Fig. |2fd)). Unlike the peaks in 
AI 2 (s) associated with crossovers in I(s), the width and 
height of these secondary peaks are r-dependent. 

Rare trajectories and dot occupation — We now return 
to the regime with /ii.2 < e and consider the additional 
peaks in AI 2 (s), which occur along the lines \i\ = sT 



and s = — u 2 /T (see Fig. [2jd)). While these features 
necessarily originate from the form of the LD function 
0(s), or the dynamical order parameter I(s), a physical 
explanation of their existence is not obvious. However, 
we gain insight by studying the state of the dot as a func- 
tion of the s-bias, rid(s), given in Eq. As shown in 
Fig. [2jf), nd(s) shows a marked crossover from zero to 
0.5 which coincides with these peaks in A/ 2 (s), showing 
that the dot occupation becomes correlated with features 
in the FCS. This crossover in dot occupation results from 
a crossover in the tunnelling rates on and off the dot as s 
and /ii.2 are tuned. In Fig. [2je) we see at the dot occu- 
pation nd(s) is s-independent when T = 0, and therefore 
does not give rise to features in A7 2 (s). 

Conclusions — We have presented a new approach to 
the study of FCS based on a thermodynamic formalism, 
taking the fully-coherent dynamics of the resonant-level 
model as an exactly-solvable example. We identified close 
analogues of equilibrium phases with sharp crossovers be- 
tween dynamical phases marked by large peaks in the 
fluctuations of transported charge. The s-biased current 
I(s) contains all cumulants in the FCS through its s- 
derivatives. However, studying FCS this way reveals a 
remarkable mapping between the FCS of quantum-dot 
systems across parameter space. We further show that 
the FCS may be mapped to the current through quan- 
tum dots at several values of the dot-lead coupling and 
applied bias. Features in the current fluctuations become 
correlated not only with the behavior of the current, but 
also with the occupation of the dot in rare trajectories, 
rid (s)- At finite bias, the dynamical phases are marked 
by sharp crossovers rather than transitions. However, we 
anticipate that transport systems with more complex in- 
ternal dynamics will show sharp dynamical transitions in 
this space of rare trajectories. We believe that analysis 
of more complex systems with the approach developed 
here, for example those with interactions [TT] or dissipa- 
tion [25], will prove fruitful. 
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SUPPLEMENTARY MATERIAL 

for 

Dynamical Phases in the Full Counting Statistics of the 
Resonant-Level Model 

In this supplementary material, we sketch the deriva- 
tions of 0(s) and n^s) for the resonant-level model 
and identify the mapping of s ^ statistics in the 
resonant-level model to the unbiased statistics of a more 
general class of models. 

Derivation of 6(s) and n<j(s) — From the definition of 
6(s) in Eq. Q, the derivation of Eq. ^ follows Bernard 
and Doyon [13] but with a rotation in the complex s- 
plane. We sketch the details as we use these as a founda- 
tion for our derivation of the dot occupation nd(s). We 
begin by expressing the dynamical partition function as 



1 ( f 2n 
Z t (s) = — Tr yp ss J drj e 



F(s,t,ri) 



where 



,F(B,t,rj) _ (s/2-M?)Q -sUlQUt Js/2+ir,)Q 



(9) 
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The dynamical partition function may be evaluated us- 
ing a formula by Klich [7] Ti{e A e B e c ) = Det(l + e a e b e c ) 
where A, B, C are operators bilinear in fermion operators 
and a, b, c are associated single-particle operators. Con- 
structing single-particle matrix elements as 2 x 2 blocks 
in the even and odd mode operators, with single-particle 
states of energy p denoted by \p) = (a£|vac), &£|vac)) T , 
one expresses the matrix elements (p\ef( s > t, ' n ' — l\p') = 
M(p, r])(e l ( p ~ p — l)/ip, neglecting exponentially decay- 
ing terms when t 3> gvpT~ 2 . Here, / is the single-particle 
operator corresponding to F and 



M(p, rj) 



-2 sinh - 
2 



T p (T p — l)(cr y cos 77 + a z sin rj) 



T, 



cosh ■ 



I sinh ■ 



(11) 



x,y,z are the Pauli matrices with I the identity. 



Here, o~ 

From the Klich formula we then have 

r2 * dr, Det[J + r 5S e^ s >*>'?)] 
'a 



Zt(s) 



2n 



Det[J- 
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dr) 
2^ 



Bet[I + N{e n>s ' t ' 7l) -1)] (12) 



where r ss is the single-particle operator corresponding to 
the Hcrshfield density matrix p ss and N = (I + r ss )~ 1 r ss 
which is diagonal in p with elements N p = (p\N\p) = 



,(-) 



a x . With careful evaluation of log Z t (s), one 
ultimately finds this to be linear in t with contributions 
only from diagonal blocks in p 13J so that Eq. ([5| is 
found. 



We now turn to the derivation of n<j(s), from the defini- 
tion in Eq. which is plotted in Fig.[2]je) and (f). Since 
d)d commutes with projection operators on the Hilbert 
space of the leads, we find 
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where our construction allows the Klich formula to be 
used. Thus 

f 271 dn d , , 
n d( s ) = I logDet 



I + Ne f( S ,t,v) e D(t) 



2n dj - J j=o 

t t (14) 

where e D ^ is the single-particle version of e jC/t <^du t _ 
Now, for brevity, denoting M(p, rf) = M p and using 
5 t {p-p') = ( e l (P-p')* - l)/i(p-p'), we find 

M p S t (p-p') + (e* - l)J pp , e ^-f')* + 

(e* - l)M p J pp/ e^-P'^S t (p - p') (15) 

where J p w = T~ 2 w p w*,(I + <J z )/2. Using this and ex- 
panding ( 14 ) we find 
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Equation (16) proves to be time independent. To see 
this, we use the identity 8 t (p) = f*dt'e ipt fr om which 
we see the k th term in the sum after ^'-differentiation has 
the form 



J dt\... dtk-i J dpi . . . dp k 



n>d(s) = J dti 

k k 

jt e ip m {t m -t m -- L ) e -i{ V i-p l+1 ){t l -t) 

1=1 m=l 

\^Y.U + e ^U bPiStfa - pi+i)] (17) 

where we consider, instead of the specific pj -dependent 
matrices, general distributions over ctj and j3j as in |13j 
and omit the ry-integration. Repeatedly using the Fourier 
representation of the delta function, we find the above ex- 



pression collapses to k J dp 



Therefore Eq. ( 16 ) takes the explicitly time-independent 
form 
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with tr denoting the trace over 2x2 blocks. After re- 
summing the series, we find 
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which, upon substituting the 2x2 matrices into the integrand and performing the 77-integral yields 



n d {s) 



dp \w p \ 2 
2^ 2 



2n 



(+) 



(4 +) - 2n p 1) n p 2) )T p (l - cosh s) + 2n p } T p sinh s - 1 



(20) 



where n 



(±) 



(1) 1 (2) 

n p ±n p . 



Mapping s ^ statistics at finite temperature to s = 
statistics — In the main text of the paper, we highlight 
how in the zero temperature limit the rare, i.e. s / 0, 
charge trajectories of one system are the typical, i.e. 
s = 0, trajectories of another. This can be seen by a 
simple rescaling of the coupling r. At non-zero temper- 
ature, the order parameter I(s) is still an integral over 
all lead modes I(s) = J I p (s)dp/2ir as in Eq. Q, but the 
integrand I p (s) is modified by s and takes the form 



We note that T p = T p when T p = or 1 and, for any real 
s, < T p < 1 implies that < f p < 1. 

In general, the transmission coefficients take a new p- 
dependent form which cannot be captured by a model 
with energy-independent couplings between the dot and 
the leads. However, discussed in the main text, an 
approximate mapping exists when \pk\i\jT 3> 1, with 
pi < < pi (e = 0), which avoids energy-dependent 
couplings. As above, the form of I(s) is an integral over 
all p-modes. However, where p 3> fJ-2 the integrand I p (s) 
is approximately 



"n^ng - 1) + e s 4 2) (l - n p 1] )]T p 



[e- s 4 1} (4 2 -l) +e s (n ( P ) -l)n p 2) + (2n ( p +) ~2n ( p l] n (2) )}T p -1 

(21) 

From this expression and its derivatives we find all the 
cumulants of the s-biased system. These are the same as 

those of an unbiased system with lead occupations 

(2) 

and n p ' and effective transmission coefficients T p . We 



4(e(p-^i)/ T + l)p 2 + (eiP-^/T+s _|_ i) r 4 



(25) 



find 



with r s = e -s / 4 r as in the main text. It is clear that 
this expression that the form of the s = integrand, but 
with r — ¥ t s and an effective modification to the chemi- 
cal potential fix in the second term in the denominator. 
An equivalent expression exists involving only the lower 
chemical potential /i 2 when p <C fJ,i- 



(1) - 
n p e 



(e- - 1) 



(22) 



and 



(2) s 

n v ■ e 



i + 4 2 V - 1) 



(23) 



which amounts to shifts in the chemical potentials in the 
two leads where p\ — > p\ + sT and /x 2 — » P2 — sT. (This 
is consistent with our understanding that the s-bias on 
the statistics has characteristics of a real potential bias; 
this is discussed in the main text with reference to the 
s-bias s/ = o required to bring I(s) to precisely zero.) The 
transmission coefficients must also be modified in a way 
that is more involved than at T = 0. The new trans- 
mission coefficients T p take the form of Eq. ^ but with 
p-dependent parameters r — > f p = t\ p with 



1 + (e" s - l)nW 1 + (e s - l)n£ 



,(2) 



(24) 



4(e-(P-M2)/T + ^ p 2 + ( e -(p-„ 2 )/T-s + ^ T 4 



(26) 



However, provided \fix t 2\/T ^ 1, when /i 2 p <C /Ui, we 
find the simple /ii^-independent form I p (s) — r^/(4p 2 + 
Tg). Therefore we are able to find a new system, with 
modified parameters r' and p[ 2 , such that I p (0) for this 
new model is almost exactly equal to I p (s) for the first 
system. We do this by matching the regimes p /ii and 
p P2 which, due to the /ii 2 -independent form for in- 
termediate p, leads to I p (0) for the new system matching 
I p (s) for all p. For example, at large p we require 



4( e (P-^i)/ T + l)p 2 + ( e (P-^)/ T + s + l)r s 4 



4( e (p-Mi)/T + l) p 2 + ( e (p-Mi)/T + 1) 7 



(27) 



from which we require the new system to have r' = r s and 
upper chemical potential given by fx[ = p\ in Eq. ^ . p! 2 
is found by an equivalent condition, with the agreement 
becoming exact in the limit T — > (see Fig. 3(b) in the 
main text). 



